**Volume Title** 

ASP Conference Series, Vol. **Volume Number** 
**Author** 

© **Copyright Year** Astronomical Society of the Pacific 



The Odd Meanderings of the IMF Across Cosmic Time 

Romeel Dave 1 

1 Astronomy Department, University of Arizona, Tucson, AZ 85750 

Abstract. It is difficult to reconcile the observed evolution of the star formation 
rate versus stellar mass (SFR-M*) relation with expectations from current hierarchical 
galaxy formation models. The observed SFR-M, relation shows a rapid rise in SFR(M„) 
from z — — » 2, and then a surprisingly lack of amplitude evolution out to z ~ 6+. 
Hierarchical models of galaxy formation match this trend qualitatively but not quanti- 
tatively, with a maximum discrepancy of ~ x3 in SFR at z ~ 2. One explanation, albeit 
radical, is that the IMF becomes modestly weighted towards massive stars out to z ~ 2, 
and then evolves back towards its present-day form by z ~ 4 or so. We observe that this 
redshift trend mimics that of the cosmic fraction of obscured star formation, perhaps 
hinting at a physical connection. Such IMF evolution would concurrently go towards 
explaining persistent discrepancies between integrated measures of star formation and 
present-day stellar mass or cosmic colors. 



1. Introduction 

There is currently much debate on whether the stellar initial mass function (IMF) is ev- 
erywhere and everywhen invariant, as highlighted by this excellent conferen ce. Direct 
evidence for IMF variations are scant and controversial (Bastian et al. 2010, and con- 



tributions to these proceedings by Meurer, Lee, and Hoversten). One possibility is that 
the cluster-intrinsic IMF is invariant, but sampling and/or truncation effects convolved 
with a varying cluster mass distribution yields a varying IMF when averaged over en- 
tire galaxies (see e.g. contributions on the IGIMF theory). Theoretically it is difficult to 
envision that, despite enormous variations in the physical cond itions of star f ormation, 
the IMF somehow conspires to remain completely invariant dKroupal [20021) . But the 
sense in which the IMF should vary is not well understood, and hence more empirical 
constraints are needed. 

Recently, a cosmological approach to constraining the IMF has gained traction. 
Here, the evolutionary properties of galaxies are studied to test whether they are con- 
sistent with an invariant IMF. Early efforts involved comparing the observed integral of 
cosmic star formation (which traces massive star evolution) with the present-day stellar 
mass density (which traces ~ IM Q stars). Such studies typically favored a cosmically- 
averaged IMF th at was weighted towards more massive stars than t he present-day Milk y 



Way disk IMF dMadau et all 1 19981: iBaldrv & Glazebrookl 120031 : IWilkins et ailb oOS). 



Other approaches include studying the color evolution of passive galaxies (Ivan Dokkuml 
2008), and applying ad ditional constraints from extragalactic background light mea- 
sures dFardal et al.ll2007l) . In general, such studies again favor an IMF more w eighted 



towa rds massive stars to satisfy all constraints. That said, a well-known theorem (Ba stian et al 



20101) states that any problem in galaxy formation can be resolved through a suitable 
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choice of IMF, and hence such arguments are often dismissed as the last resort of a 
scound rel. This fail s to deter us theorists, who have been called far worse. 

In lDavel(l2008h . we argued that the evolution of the star formation rate-stellar mass 
(SFR-M*) relation may be suggestive of a time-varying IMF for typical star-forming 
galaxies. At that time, reliable observations of SFR- M* only extended out to z ~ 2. 



Improved data has now co nstrained SFR-M* to z ~ 7 (IStark et all 120091 : lLabbe et al. 



20101 : Gonzalez et al.ll2010h . providing an opporunity to extend such arguments to very 



high redshifts. Concurrently, Herchel has provided improved bolometric constraints on 
star formation in galaxies out to z ~ 2, mitigating some of the systematic uncertainties 
that may have affected earlier works. Theoretical work has als o progressed, as to y 
models have been developed that can explain SFR-M* evolution ( Bouche et al.ll20l6T) . 
at the cost of introducing arbitrary modifications to galaxy star formation histories. 
With such fervent recent activity, the time is right to revisit our previous arguments. 
This is the goal of these proceedings. The end result, we argue, is that the evidence still 
favors some discrepancy in the SFR vs. M* evolution in galaxies that could be resolved 
by a time-varying IMF, though to do so requires a modestly varying IMF that evolves 
non-monotonically with redshift. 



2. The SFR-M* Relation 



The SFR-M* relation has recently received much attention as an important barometer 
of galaxy evolution theory. Models, particularly cosmological hydrodynamic simula- 
tions, have long predicted that there should be a tight and slightly sub-linear relation 
between these quantities (e.g. iDave et al.ll2000l: jFinlator et al.ll2006l : iDave et al.ll2006l) 
that evolves slowly upwards with redshift (iDavel 20081). Observations over the last fiv e 
years h ave confirmed thes e trends out to z ~ 1 dNoeske et al]|2007l : lElbaz et ai1l2007l) . 
z ~ 2 dDaddi et all 120071) . and beyond (Gonzalez et al.ll2O10h . In models, this trend 
arises because cold and smo oth accretion is the dominant accretion path for fuelin g 
star formation at all epochs dMurali et al.ll2002l : iKeres et al.l 120051 : iDekel etaii r2009): 
the tightness of SF R-M* empirically constrains the contribution from starbursts (e.g. 
Noesk e et al. I l2007h . The broad agreement between theory and observations is a key 
success of current galaxy formation models. 

In detail, however, it has been noted that SFR-M* predicted in models (e ither hy- 
drodynamic or semi-ana l ytic) fails t o match the observed amplitude evolution (Elbaz et al 
20071 : iDaddi et alJl2007l : lDavell20 08). This is shown in Figured] which shows the evo- 



lution of the SFR-M* relation at a fiducial stellar mass of M* 
cosm ological hydrodynamic simulations (details of which are in lOppenheimer et al. 
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2010h . The green curve shows results from a simulation including outflows with scal- 
ings as expected for momentum-driven winds; this outflow model has proved successful 
at matching a wide range of galaxy and IGM data from z ~ - 6. The blue curve shows 
a simulation with no galactic outflows included , which strongly overpredicts the global 
stellar mass formed dOppenheimer et aLll2010l) . but nevertheless produces a SFR-M* 
relation that is not markedly different from the strong outflow case. This illustrates 
that despite dramatic variations in predictions owing t o different feedback mode ls, the 
SFR-M* relation is relatively insensitive to feedback dDavell2008l : lDuttonll2010t) . The 
basic reason is that feedback suppresses star formation which likewise suppresses stel- 
lar mass growth, moving galaxies downwards along the (nearly linear) SFR-M* rela- 
tion. 
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Figure 1. 
Z = - 7 



The amplitude evolution of the SFR- M, relation at Af» = 10 I0 M n from 
Observations (circl es) are shown from lElbaz et alj j2007h at z = 0, 1, 
and from iGonzalez et all (1201 Oh at z > 2. Dashed lines indicate +0.3 dex spread 
around the mean relation. Red dot roughly indicates the downward connection of 
UV-derived SFRs based on Herschel data. Simulation results are shown for our 
favored run employing momentum-driven winds (green solid) and a run with no 
galactic outflows (blue dotted). 



Observations from lElbaz eTall (|2007l. atz = 0, 1) and lGonzalez et al.l(|2010l. z > 2) 
are shown as circles, and a typical lo~ observed scatter around the relation of ±0.3 dex 
is shown by the dashed lines. Note that this is not the statistical or systematic uncer- 
tainties in M* or SFR, which are typically smaller than 0.3 dex at low-z but could be 
significantly larger at high-z. As can be seen, the observed SFR-M, relation amplitude 
matches well at z ~ 0, and then the observed SFRs are increasingly larger going out to 
z ~ 2. At this point, there is an abrupt change in the trend, and the observed amplitude 
is essentially invariant out to z ~ 7. The simulations broadly predict a similar trend, but 
do not reproduce the abrupt change at z ~ 2. 

The discrepancies are modest, especially considering the uncertainties involved, 
maximizing at ~ x3 at z ~ 2. But what is particularly puzzling is that at z ~ 2, the 
observed star formation rates can exceed the accretion rate into galaxy halos, which 
is difficult to accomodate in current models that favor continual gas supply ("cold 
streams") to fuel star formation. The implication is that galaxies must store up large 
reservoirs of gas, and then something triggers its rapid conversion into stars. If star for- 
mation proceeds exponentially, then galaxies quickly move upwards along the SFR-M* 
relation iMaraston et all (120101) . preserving the (nearly) linear slope at a fixed ampli- 
tude. Current hierar chical galaxy formation models do not yield this behavior naturally. 
Bouch e et al. (2010) noted that an ad hoc way to accomplish this is to explicitly prevent 



all star formation in halos with M < 10 n M o , which accumulates a large reservoir that 
is then rapidly consumed once the threshold halo mass is crossed. But this is difficult 
to arrange physically, and the lack of any discernible transition in the galaxy popula- 
tion around these halo masses (as compared to, say, around 10 12 M o ) argues against this 
scenario. Hence matching the z ~ 2+ SFR-M* relation appears to require a fundamen- 
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tal modification to a galaxy formation paradi gm that is oth erwise quite successful at 
reproducing a wide range of observations (e.g. lBensonll201fj|) . 

An alternative possibility is that there ar e significant systematics in observational 
measures of SFR and/or M*. As discussed in [Pavel (120081) . it is easier to accommodate 
overestimated SFR's than underestimated M*, since stellar mass is cum ulative and there 
are lower-z constraints that must be satisfied. The iDaddi et all (120071) z ~ 2 data rely 
on extinction-corrected UV SFRs, which recent Herschel data suggest overestimate the 
true (bolometric) SFRs measured from far-IR data by 50% or perhaps up to ~ x2. 
Such a reduction is indicated by the red point on Figure Q] This would alleviate the 
discrepancy somewhat but not completely, and in a sense makes it harder to reconcile 
since a key systematic has been mitigated. We note that 24//-based SFRs have been 
shown to be very poor at z ~ 2 sinc e PAH emiss ion that dominates the rest-frame Sfi 
emission is poorly understood (e.g. lKelsonll2010l) . but such SFRs are not the basis for 
the SFR-M* observations shown. Nevertheless, better data is required to pin all the 
systematics down, and a proper accounting may yet remove the discrepancy altogether. 

Finally, Figure Q] shows that feedback does have a noticeable impact on the shape 
of SFR-M* evolution. In particular, the momentum-driven wind scalings suppress 
star formation in smaller halos more, resulting in less early star formation in better 
agreement with data. One might envision that an even stronger feedbac k model could 
compl etely reconcile the theory and observations. However, the model of iBouche et al. 



( 20101) indicates that the feedback required is quite extreme, essentially suppressing 
all star formation even in fairly sizeable halos. Our experience is that models with 
significantly stronger feedback have substantial difficulties producing enough early star 
formation to match data on z ~ 6+ galaxies, the mass-metallicity relation, and the 
metallicity of the IGM at early epochs. Perhaps some compromise will eventually 
be found between theory and observations, and certainly the downward revision of 
observed z ~ 2 SFRs goes towards this. This probably remains the most likely solution 
to the SFR-M* discrepancy, but it is not the most interesting one. 



3. A Time- Varying IMF? 

A more interesting and far-reaching systematic that could lower inferred SFRs is if the 
IMF is weighted towards more massive stars at high redshift, either by being top-heavy 
or bottom-light. Here I use top-heavy to indicate a changing upper-end slope, and 
bottom-light to refer to the suppression of low-mass stars. If this is to be the solution, 
then the IMF must evolve with redshift towards more top-heavy/bottom-light at z ~ 2, 
and then reverse its trend back towards a normal IMF by z ~ 4+. At z > 5 the systematic 
effects in M* and SFR are currently sufficiently large that it would be premature to infer 
anything from the small discrepancies at those redshifts. 



One possible form of IMF variation, highlighted in Dave (2008) and Ivan Dokkum 



(120081) . is that the IMF characteristic mass evolves with redshift. For a K roupa IMF , 
Le - l^M x M '° 3 for M < ^imf and oc M-° for M > Mjmf, [DayJ (|2008|) 

showed that one evolutionary form that reconciles SFR-M* evolution out to z ~ 2 is 
Mimf = 0.5(1 + z) 2 M Q . However, beyond z ~ 2 the IMF must evolve back towards the 
present-day form, and moreover Herschel data has shifted the z ~ 2 point downwards. 
A general form that would accommodate these trends is 

M IMF = 0.5(1 +z) 3 -° J5z M Q (1) 
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Figure 2. The IMF characteristic m ass Mmv in various IMF evolutionary forms. 
The solid black line shows the form of Dave (2008), which reconciles SFR-M, from 
Z = — > 2. Recent results have suggested the IMF returns to its present-day form 
by z ~ 4, motivating a different evolutionary form. The green dashed line, with 
Mimf = 0.5(1 + z) 3 ~ l)J5z M B , follows the old form out to z ~ 1.5, and then returns 
towards the present-day IMF by z — 4, broadly accommodating current constraints. 



from z. = — > 4, remaining at 0.5M o at z > 4. As shown in Figure P21 this roughly 
mimics the Mimf evolution of lDava (120081) out to z ~ 1-5, and then turns over, reducing 
the characteristic mass by almost a factor of two at z ~ 2, and eventually returns to 
the present-day value at z - 4. This is purely an eyeball estimate, and furthermore 
this choice of redshift dependence is merely a convenient parameterization; there is no 
physics in it. In reality, the IMF should be tied to some physical property of galaxies 
that evolves with redshift, such as star formation rate or surface density. 

The SFR-M* relation does not constrain the detailed form of IMF variation. Other 
plausible ways to reconcile the SFR -M* would be to make the upper-end slope of the 
IMF more shallow with time (e.g. IWilkins et al.ll2008T) . In detailed studies of lensed 
Lyman Break (z, ~ 2 -3) galaxies, the u pper end of the IMF does not appear to be 
substantially different (Oui der et al.| [2009) though current constraints may not rule out 
mild changes that may resolve the SFR-M* discrepancy. Also, constr aints on t h e IMF 
in z ~ 2 sub-millimeter galaxies indicate that it is consistent with a Chabriei] d2003h 
IMF, perhaps slightly favoring a more top-heavy one, but clearly inconsistent with dra- 
m atic departures from a Salpeter slope as postulated in e.g. the semi-analytic models 
of iBaugh et al.1 (120051) . On the other hand, the ratio of Ha to UV emission in local 
star- forming galaxies co r relates strongly with its star formation rate and/or surface den- 
sity ( Meurer et al.ll2009l : iLee et al.ll2009h . for which one interpretation is that the IMF 
has a shallower upper-end slope where there is more concentrated star formation. Since 
z. ~ 2-3 galaxies ubiquitously show more concentrated star formation, this could empir- 
ically imply a top-heavy IMF in such systems. It remains to be seen if simply applying 
constraints from local galaxies would yield the required evolution owing purely to the 
size and star formation rate evolution of typical galaxies, but this should be investigated. 
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An IMF that returns to normal at very high redshifts is also preferred by avail- 
able constraints at those epochs. For instance, current models assuming a normal IMF 
are able to form en ough stellar mass in very early (z ~ 6+) objects as observed (e.g. 
Finl ator et al. but a significantly bottom-light or top-heavy IMF would create dif- 



ficulties. Another early-universe constraint comes from old globular clusters formed at 
high-z, which are consistent with a present-day IMF. Hence there is a growing concen- 
sus that the IMF in the early universe is actually quite similar to today's. Of course, Pop- 
ulation III (i.e. very metal poor) stars are expected to have a quite top-heavy IMF, but 
it is unlikely tha t these contribute sig nificantly to early galaxies observed to date (e.g. 



Dave et al. 2006; Finlator 



ntnbut e sig n 
et alJl201Ch . 



The central question is, what could cause the IMF to depart away from the present- 
day IMF in such a fashion? The non-monotonic evolution is particularly challenging to 
understand. This is a more a question for star formation theorists, but there is one pos- 
sible intriguing connection: The required IMF evolut ion roughly mimics th e cosmic 
fraction of obscured star formation as determined by [Bouwens et al.1 ([2009, see their 
Figure 11), which peaks at z ~ 2 - 3 and falls off to higher and lower redshifts. Hence 
if there was something about an obscured mode of star formation that systematically 
caused the IMF to be weighted towards massive stars, this may yield the desired trend. 
A possible local test would then be to look in nea rby ULIRGs, where som e controver- 
sial claims of a top-heavy IMF have been made (Rie ke et al]|l980l Il993l) . Obviously 
this connection is highly tenuous, but may provide an avenue for further exploration. 



4. Cosmically-Averaged Star Formation 

The IMF variation required to reconcile theory and observations of SFR-M* also broadly 



reconciles the observed c osmic star formation and stellar mass growth histories (Dave 
2008 ; IWilkins et al.ll2008l ). and extragalactic background light constraints (IFardal et al. 
2003). Clearly these quantities are all interrelated, so it isn't terribly surprising that 



the same IMF fixes all of them, but it provides a useful consistency check. Another 
approach (see Wilkins, these proceedings) is to use the cosmic SED to constrain the 
IMF; once again, this shows a mismatch between the integrated cosmic star formation 
history and the present-day mean cosmic color in the analogous sense that there is too 
little near-IR light relative to the UV/optical light together with expectations from the 
observed cosmic star formation history. Present-day integrated measures do not neces- 
sarily require an evolving IMF, but rather only a cosmically-averaged IMF that differs 
from the local one. However, it would be puzzling if our locally-measured IMF was 
different than the (invariant) IMF over the rest of cosmic time. 

Conversely, differential measures of stellar mass growth and star f ormation rate 
could imply an evo lving IMF, and some claims to that end have been made (Wil kins et al.l 
20081 : lDavell2008h . Howeve r, such claims are cont roversial. A particularly strong ar- 



gument was forwarded by iReddy & Steidell (120091 also see these proceedings), who 



noted that by making careful corrections for dust and incompleteness, it is possible to 
reconcile the global star formation rate and stellar mass density evolution at z ~ 2. 
This would not counter arguments based on individual galaxies such as the SFR-M* 
discrepancy, but would make observations of cosmic star formation and mass growth 
internally self-consistent at z ~ 2. 

In actuality, how ever, this argumen t ends up merely pushing the problem to lower 



redshifts. To see this, Red dy & Steid el (2009) notes that their incompleteness correc- 
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tions result in 57% of the present-day stellar mass already having been formed by z ~ 2. 
Conversely, current observations of the cosmic star formation hist ory imply that only 
~ 30% of total cosmic star formation occurs prior to z = 2 (e.g. iHopkins & Beacom 
2006; Far dal et aT1l2007 ) — this is independent of assumed IMF, so long as it is time- 
invariant. Hence the growth from z = 2 — > of stellar mass does not match expectations 
from the observed star formation history, in the sense (once again) that star formation 
rates are measured to be too high. Hence one can "fix" the problem at z ~ 2 and have 
difficulty with evolution to z = 0, or one can "fix" the problem at z = and then have 
difficulty reconciling z ~ 2 observations. This argues for systematics in SFR or M t 
measures that change with redshift, be they related to the IMF or not. 



5. Conclusion 

There is still no definitive evidence for IMF variations across cosmic time. That said, 
the issues motivating serious consideration of such variations have yet to be fully re- 
solved. New data on the SFR-M* relation out to z ~ 7 suggests that, if IMF varia- 
tions are to explain the discrepancy with theoretical expectations, deviations from the 
present-day IMF must maximize at z ~ 2 and then lessen to higher redshifts, returning 
to a standard IMF by z ~ 4-5. A possible form that would roughly yield this is a Kroupa 
IMF with a characteristic (turnover) mass evolving as M 1MF = 0.5(1 + z ) 3 ' OJ5z M Q . 
While the required non-monotonic IMF evolution is not straightforward to understand 
physically, a quantity showing analogous behavior is the cosmic fraction of obscured 
star formation, hinting at a possible physical connection. 

The form of the IMF variation remains difficult to pin down. The best (albeit 
controversial) evidence for local variation is in the form of a changing upper-end slope 
with star formation rate. The high star formation rates of high-z galaxies would then 
naturally result in a more top-heavy IMF. But direct observations of z ~ 2 - 3 galaxies 
disfavor significant changes in th e upper-end IMF s lope. Meanwhile, t he evolutionary 
form suggested by iDavel ([2008) (or, similarly, by Ivan Dokkuml 120081) towards more 
bottom-light with redshift is more difficult to test directly. It is also possible that the 
intrinsic IMF is invariant, but the larger typical sizes of molecular clouds in high-z 
galaxies results in a galactic-averaged IMF that yields proportionally more massive 
stars. An IMF weighted towards massive stars at intermediate redshifts would also go 
towards resolving the persistent discrepancies between integrated measures of cosmic 
star formation history and present-day measures of stellar mass or cosmic light. 

In the near future, surveys such as the Spitzer Extragalactic Deep Survey (PI G. 
Fazio) and the Hubble CANDELS survey (Pis S. Faber and H. Ferguson) will provide 
unprecedented multiwavelength samples of galaxies out to z ~ 6 and beyond. Although 
this will not necessarily pin down systematic effects in SFR or M* measures, it will 
enable tests of different galaxy formation scenarios with greater statistical precision. 
For instance, the "reservoir" scenario where gas consumption proceeds exponentially in 
galaxies predicts a significantly different evolution to high-z of the stellar mass function 
and halo occupancy than hierarchical galaxy formation models. Hence these surveys 
will be critical for testing current galaxy formation models, thereby indirectly testing 
ideas such as IMF evolution. We look forward to engaging in such comparisons. 
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